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SUMMARY 


A numerical study of computations in backward- facing steps with flow 
separation and reattachment, using the Reynolds stress closure is 
presented. 

The highlight of this study is the improvement of the Reynolds - 
stress model (RSM) by modifying the diffusive transport of the 
Reynolds stresses through the formulation, solution and subsequent 
incorporation of the transport equations of the third moments, u ^ u j u ^ c » 
into the turbulence model. The diffusive transport of the Reynolds 
stresses, represented by the gradients of the third moments, attains 
greater significance in recirculating flows. The evaluation of the 
third moments by existing algebraic correlations is inadequate but the 
third moments evaluated by the development and solution of the complete 
transport equations are superior than those obtained by the algebraic 
correlations . 


A low-Reynolds number model for the transport equations of the third 
moments is developed and considerable improvement in the near-wall 
profiles of the third moments is observed. The values of the empirical 
constants utilized in the development of the model are recommended. 

The Reynolds -stress closure is consolidated by incorporating the 
equations of k and £, containing the modified diffusion coefficients, 


and the transport equations of the third moments, u.u.ul, into the 

1 j K 

Reynolds stress equations. Computational results obtained by the 


original k-c model, the original RSM and the consolidated and modified 
RSM are compared with experimental data. Overall improvement in the 


i 



predictions is seen by consolidation of the RSM and a marked improvement 
in the profiles of is obtained around the reattachment region. 
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CHAPTER 1 


INTRODUCTION 


1.1 Characteristics of Backward-Facing Step Flows 

The phenomena of flow separation with the consequent wake formation, 
reattachment and recirculation occur in a vast array of engineering 
applications. Engineers use it to their advantage in augmenting the 
performance of equipment like combustors, heat exchangers and reactors. 
On the other hand, they also devote considerable effort in streamlining 
airfoils and similar boat-tailed bodies in order to minimize the 
incidence of flow separation. Very often this phenomena may be held 
responsible for the optimal operation of a system, which just goes to 
show the enormous significance of being able to correctly study and 
model it. 

Two-dimensional backward- facing step flows have been extensively 
used in experimental as well as theoretical (numerical) methods to study 
the phenomena of flow separation, recirculation and reattachment. This 
class of flows possesses all the important aspects of separation, 
recirculation and reattachment, but at the same time, owing to their 
geometrical simplicity, they render themselves to easier experimentation 
and computation, without loss of flow complexity or susceptibility to 
modeling modifications. 

Moreover, due to greater interest in this field, there exist 
relatively large sets of data on such class of flows which facilitates 
easier validation with modeling trends. However, some of the data like 
measurements of higher order moments of turbulence have not yet been 
explored widely. 


Figure 1.1 depicts the flow field and the regimes in a typical 
backward- facing step. The upstream boundary layer separates at the 
sharp corner forming a free -shear layer. The separation line is 
straight and fixed at the edge of the step and there is only one 
separated zone. The streamlines are nearly parallel to the wall at the 
separation point, and significant upstream influence occurs only down- 
stream of separation. 

An ordinary boundary layer develops on the upper wall separated by 
an inviscid core from the viscous and separation effects on the lower 
wall. In case of an open top geometry, however, the inviscid core 
(potential flow) extends upto the top forming a free shear boundary at 
the surface. 

The separated shear layer appears to be similar to an ordinary 
plane -mixing layer through the first half of the separated flow region. 
The dividing streamline is slightly curved and the shear layer is not 
influenced by the presence of the walls. The important difference 
between the separated shear layer and the plane -mixing layer lies in the 
fact that the flow on the low speed side of the shear layer is highly 
turbulent due to the flow reversal and recirculation, unlike the low 
turbulence levels encountered in plane-mixing layers. 

The separated shear layer curves downward in the reattachment region 
and impinges on the wall. A fraction of this fluid gets deflected 
upstream into the recirculation region due to strong adverse pressure 
gradients causing the formation of a chaotic primary two-dimensional and 
a small three-dimensional eddy inside this region. Armaly et al. 6 
report the existence of two recirculation regions on the step-side wall 
(inside and outside this recirculation zone), and one on the upper wall. 



The length of the recirculation region varies 6,32 ’ 56 periodically; 
this is attributed to a low frequency "flapping" motion of the shear 
layer but in the experiment of Ruderich and Fernholz* 5 , however, no such 
flapping was observed. 

The shear layer is subjected to the effects of stabilizing 
curvature, adverse pressure gradients and strong interaction with the 
wall in the reattachment region. The large eddies are confined by the 
wall which also results in their distortion by effectively irrotational 
mechanisms resulting in a spectacular decrease in the Reynolds stresses 
and the third moments thereby leading to the transport of turbulent 
kinetic energy and shear stresses toward the solid wall from the region 
of maximum turbulent intensity. 

After impinging on the wall, the remnant fluid in the reattached 
shear layer continues to flow downstream developing a boundary layer on 
the wall resulting in the formation of the redeveloping layer. 

Downstream of reattachment, the Reynolds stresses and the third 
moments continue to decay rapidly for a distance of several step 
heights. Measurements of Bradshaw and Wong and Smyth reveal that the 
outer part of the reattached shear layer still has most of the 
characteristics of a free shear layer as many as 50 step heights 
downstream of reattachment verifying the persistence of large scale 
eddies which are developed in the separated free -shear layer. 

The process of separation and reattachment in backward- facing step 
flows, hence, makes them rather complex, and thus much different from 
boundary layer or mixing layer types of flows. Therefore, the 
turbulence model used for predicting this class of flows should not just 
be capable of accounting for all these phenomena but should also possess 
the sophistication of being accurate in doing so. This led to the 


modification and development of the Reynolds - Stress Model (RSM) by 
incorporating the third moments into it as shown in the subsequent 
sections . 

1 . 2 Literature Survey 

Turbulence models capable of predicting complex flows like those 
occuring inside a backward- facing step flow have been sought by many 
researchers. Among the existing models of turbulence, the three most 
popular ones are (i) k-e Model, (ii) Algebraic-Stress Model (ASM) and 
(iii) Reynolds -Stress Model (RSM). 

These models have been used to predict 

separating and recirculating flows to an acceptable degree of 
satisfaction, but owing to the simplifying assumptions made in their 
development, these models tend to be approximations of the real physical 
phenomenon, predicting results which may not altogether be too 
realistic . 

30 

The k-e Model of Jones and Launder relies on the assumption of an 

isotropic turbulent viscosity in the solution of the momentum equations 

and the equations of k and e . As pointed out by Driver and 
19 

Seegmiller , this model overestimates the isotropic turbulent viscosity 
in recirculating flows resulting in a higher spreading rate for the 
shear layer thereby yielding a premature reattachment. 

A similar view is shared by Thompson and Whitelaw 55 who concluded 
that the instantaneous reversed flow makes a positive contribution to 
the turbulence shear stress and therefore the isotropic (eddy) viscosity 
assumption may be inappropriate in regions with flow reversals. Amano 
and Goel 5 compared the results obtained by the above three models for 
flow into a channel with a backward -facing step and found the k-e model 


to consistently underpredict the levels of the Reynolds stresses. Heat 
transfer results for flow into a pipe with a sudden expansion showed 

that the RSM improves the prediction over the k-e Model. 

A3 

The ASM of Rodi is a slight improvement over the previous model 
because it takes into consideration the presence of Reynolds stresses in 
the flow field. However, the ASM does not take into account the 
convective and diffusive transport of the Reynolds stresses but instead, 
relates it to the turbulent kinetic energy leading to algebraic 
equations for the Reynolds stresses. Modified version of this model 
containing a new dissipation rate equation wherein the production rate 
is made more sensitive to the strearawise curvature effects was developed 
by Sindir 50 . The ASM, generally, as shown by Driver and Seegmiller , 19 
results in the inability of the model to locate the peak values of the 
turbulent Reynolds stresses due to an inadequate treatment of either the 
convection or diffusion mechanism in the kinetic energy transport 
equation. 

The RSM comprises of individual transport equations for all the 


Reynolds stresses, ruu. , thus accounting for all the processes of 
convection, diffusion, production and dissipation of the stresses within 
the flow field. The diffusive transport of the Reynolds stresses is 
governed by the gradient of the third moments, u.u.u, . In order to 

i j K 

develop an expression for this, a transport equation for the third 
moments was formulated 15 and through simplifying assumptions 17,25 
various terms in it were dropped out to yield algebraic expressions for 
the third moments, in terms of the Reynolds stresses and their spatial 
gradients . 

16 

As shown by Cormack et al. , these algebraic expressions give 


satisfactory predictions of the third moments when compared with the 
experimental data for flows in channels, pipes, wall- jets and mixing 



, 27 , 36 , 29,58 

layers . 

These algebraic expressions for the third moments seemed to work 

satisfactorily for a simple boundary layer type of flows. However, as 

12 

pointed out by Chandrsuda and Bradshaw , for flow in a backward- facing 


step, the u £ u j u k Profiles beyond x/H *=12.0 had the same shape as those 
in a boundary layer but their magnitudes were larger by a factor of 10 
or more, for the same free shear velocity, indicating a need for some 
correction. They also pointed out that although streamwise decrease of 
third moments in the outer part of the shear layer roughly follows the 
decrease in the shear stress or shear stress gradient, the decrease of 
the third moments near the wall is far more rapid than the decrease of 
shear- stress or intensity gradient with streamwise distance thus 
throwing doubt upon the usefulness of the way of expressing in terms of 
a gradient diffusivity even as a pragmatic means of correlating 
experimental data. They argued that since turbulent transport by third 
moments is a very significant part of the Reynolds -stress balance in 
mixing layers, it should be modeled with reasonable accuracy in 
calculation methods implying that the algebraic relations between the 
third moments and the Reynolds-str esses (or their gradients) are likely 
to be inadequate for strongly perturbed flows and that the effects of 
approach to a wall imply that any calculation method, intended to deal 
with reattaching and separating flows, should include a fairly 
sophisticated model for the third moments, preferably based on the 
transport equations of the third moments. 

Parallel with the development of theoretical models, extensive 
experimental observations were being made on such flows using a variety 
of experimental techniques like laser-doppler and pulsed-wire 


anemometers and hot-wire probes. 


6 , 12 , 19 , 22 , 32 , 45 , 51,52 


Nevertheless , 


attempts at modeling the transport equations of the third moments and 
solving them iteratively with the transport equations of the Reynolds 
stresses have not been encountered in the literature to date. 

1.3 Objectives of the Investigation 

The literature survey points out some of the shortcomings of the 
turbulence models and recommends different means for improvement, some 
of which are considered here. 

The objective of this investigation is to attempt to overcome some 
of these shortcomings and deficiencies and propose an improved version 
of the Reynolds -Stress Model of turbulence. In order to carry out and 
complete this study, the general features that have been aimed at and 
accomplished are mentioned below. The flow geometry employed is the 
backward- facing step. 

1. Formulation and development of transport equations for the 
third moments and the recommendation of different empirical 
constants used in the modeling process. 

2. Development of a low-Reynolds number model for the transport 
equations of the third moments in order to account for the 
viscous effects near the wall. 

3. Consolidation of the RSM, wherein, the two momentum equations 
along with the equations of k, e and the Reynolds stresses, 
containing the third moments, and the transport equations of 
the third moments are solved together iteratively. 

Converged values of the third moments obtained by the solution of 
their transport equations would be used to evaluate the diffusion rates 
of the Reynolds stresses in their transport equations, solution of which 



would enable the evaluation of diffusion rates of k and e, thereby 
reducing the dependence on an isotropic assumption in the model. 


CHAPTER 2 


TURBULENCE MODEL WITH GOVERNING EQUATIONS 


2 • 1 The Reynolds Stress Model 

In order to account for the presence of the Reynolds stresses in the 
turbulent momentum equation, it becomes necessary to model and evaluate 
these Reynolds stresses suitably. 

One of the most commonly used correlation is that of Boussinesq 
approximation wherein 


flU. 5U. 2 

- pu u - p (— i + — L) - - s pk (2.1.1) 

J c a Xj ax i 3 1J 

where is known as the "turbulent coefficient of viscosity" or "eddy 
viscosity" which is analogous to the concept of laminar viscosity, /i. 
Substitution of equation (2.1.1) into momentum equation yields 

au. a i ap i a au. au. 

— + OMJ.) - - - + [(p + p ) (— 3 1 + — 1 -)] 

at ax. J p ax. p ax. ax. ax. 

J 1 j j i 

( 2 . 1 . 2 ) 

Equation (2.1.2) is also known as the Boussinesq Effective Viscosity 
Model since (p + p^_) together constitute the effective laminar and 
turbulent viscosities. 

The eddy viscosity p^ has been further expressed as 

2 

C pk 

where is a constant (- 0.09). 
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(2.1.3) 


Since the concept of eddy viscosity is based upon the assumption of 
isotropic turbulence it is commonly referred to as "isotropic turbulent 


viscosity" . 

The Reynolds stresses can also be evaluated by modeling and solving 
the transport equation for u ^ u j explicitly. 

The transport equation for the Reynolds stresses is given as 
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(2.1.4) 


The left side of equation (2.1.4) represents the convective 
transport of the Reynolds stresses and the terms on the right represent 
Term I ■* Production 
Term II « Dissipation 

Term III ■ Pressure-Strain Redistribution 


Term IV 


Diffusion 



Generally the second term of the dissipation and the last two terms of 
the molecular transport in the diffusion term are negligible. Thus for 
a fluid of uniform density p and viscosity p, the transport equation for 
the Reynolds stresses may be expressed in the form 


D Vi 


™ P . . ■ e . . + . . + D . . 

Dt ^ ^ ^ "J 


(2.1.5) 


where 


3U . au. 

hj - -< U A ~r + ”i"k r 1 ) 


a \ Sx k 





du. 

a u. 

- 2v 

l 

l 


5x k 

a \ 

P 

au 

au. 

- - 

(— 

+ — 

P 

ax. 

j 

ax i 


( 2 . 1 . 6 ) 


u.p u.p 


au.u. 


- - t u i u i\ + -ik ' v ik 

J ax^ 1J ^ p Jk p a^ 


Closures for the terms in equation (2.1.6) have been proposed by 
different researchers as shown below. 

The production term P needs no further modeling since it consists 
of known and available components, essentially the Reynolds stresses and 
the mean strain rates. 

The closure for the dissipation term has been adopted from 
Rotta and which is also supported by several researchers, where it is 


assumed that 



12 


where e - dissipation rate of turbulent kinetic energy. 

The pressure -strain closure <f > ^ has been developed independently by 
Naot et al. 38 using a two point correlation and by Launder et al. 3 * 
where a fourth order tensor is employed to simulate the pressure- strain 
effects . 

Modified versions of the pressure-strain closure containing effects 

due to temperature and gravitational field and resistance to large 

23 35 57 

anisotropy have been proposed by different researchers ’ ' . 

Preliminary studies have revealed that the closure proposed by 
Launder et al. 3 * gives better agreement with the experimental data and 
hence their closure has been retained in this study to evaluate the 
pressure-strain effects. 

According to Launder et al., the closure for the pressure- strain 
effect is given as 


*u - - c i e ( p • ; V • V p ij • ; v> 


au. au. 2 
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where 

B 1 “ (C 2 + 8)/11, B 2 " (30C 2 * 2)/55 ’ B 3 “ <8C 2 ' 2)/11 


C 1 - 1.5, C 2 - 0.4 (2.1.9) 
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G - - u.u. 
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In order to model and approximate the diffusive mechanism of u ^ u j . 
the relative magnitudes of the various terms comprising the diffusive 
terms were investigated. 

The second and third parts in the bracket of D„ in equation (2.1.6) 
have been experimentally known to be much smaller than the third moments 


(u^u^u^) and are consequently neglected. This is seen from the 

29 

experimental data of Irvin for self -preserving jets in an adverse 
pressure gradient. The profile for the total diffusion of the Reynolds 
stresses was nearly of the same magnitude as that obtained by the term 


of the third moments (u.u.u, ), indicating that the diffusive effects 

l j K 

exercised by the pressure flutuations and viscous action are neglible 


compared with those due to u^u^u^. However, we still keep the viscous 
term in order to account for the wall effect. 


17 


Therefore , following Daly and Harlow , as explained in Section 2.5, 
the diffusive transport term is cast as 
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(2.1.13) 


The turbulence model described above has been known to give 

3 

satisfactory results for sudden expansions in circular pipes, and for 
plane backward facing step channels 5 . 

In order to get as close to the real phenomenon as possible, the 
diffusive mechanism of the Reynolds stresses should be evaluated in 


terms of the third moments, u.u.u, , explicitly. This modification and 

1 j x 


14 


improvement is discussed in the subsequent section. 

Upon comparison of results obtained by the above formulations 
((2.1.5) through (2.1.13)) with the experimental data for homogeneous 
free shear layers and near-wall turbulence, Launder et al. 3A found that 
the near-wall streamwise stress component is appreciably larger and the 
transverse component much smaller than in the homogeneous free shear 
layer. The levels of the shear stress component uv were also much 
smaller in the near-wall flow. Thus, in order to compensate for this 

34 

phenomenon, Launder et al. proposed an additional wall term, to 

be added to equation (2.1.5). The modeling of this term is similar to 
the modeling of the pressure -strain term <f> of equations (2.1.8) such 
that its presence tends to increase the anisotropy of the normal 
stresses but tends to diminish the shear stress. 

According to Launder et al., 
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where Q„ and P„ are given by equations (2.1.10) and (2.1.11) respec- 
tively and z is the net distance from the walls described as 


+ I 

z x y 


With the incorporation of equation (2.1.14) into equation (2.1.5), 
the Reynolds -Stress Equation can be written as 
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2.2 k-e Model and Its Modifications 


To determine the transport equation for k, the unaveraged transport 
equation for a single fluctuating velocity component u^ is multiplied by 
u^ and then time averaged. Since k - u ~^/2, the above operation yields 
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The term 
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is the generation rate of the kinetic energy and upon invoking the 
Boussinesq's approximation, can be written as 
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The second term on the right side of equation (2.2.2) represents the 
diffusion rate and is modeled as 
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Approximating the dissipation rate e by 


and neglecting other terms, the final transport equation for k reduces 
to 

Dk 13 fj, 3k 
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where er k represents the turbulent Prandtl number for turbulent kinetic 

energy. Likewise, the transport equation for dissipation rate e is 

obtained by differentiating the unaveraged transport equation for u. 

1 3u. 

with respect to x. and then multiplying this equation through by 2u — — 
and finally time averaging the resulting equation. The final form of 
this equation is, 

De 1 8 n de e 

[(M + — ) ] + C - G - C — (2.2.6) 

Dt p 5x. a 5x. k k 

J « J 

where and are constants equal to 1.44 and 1.92, respectively, 
and is the turbulent Prandtl number for dissipation rate. 

In the modeling of the turbulent diffusive transport mechanism of 
kinetic energy k, the third moments, 

2 

U i "I 

2 

in equation (2.2.2) are decomposed into the products of Reynolds 
stresses using Daly and Harlow's 17 approximation of equation (2.1.13). 

A further assumption of isotropy led to the simplified diffusion rates 
as shown in equation (2.2.4). 

Likewise, the turbulent diffusive transport mechanism of the dissi- 


25 

pation rate e was decomposed by Hanjalic and Launder as being 
proportional to 


de 


- u.u. 

e 1 J ax. 

J 


and a further assumption of isotropy led to the final form of equation 

( 2 . 2 . 6 ). 
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A low-Reynolds number model was developed by Jones and Launder * 
to predict laminar ization phenomenon, in which additional terms bearing 
functional viscosities were incorporated in the k- and e- equations. 

Both the high and low-Reynolds number versions of k-e Model were 
applied by Amano and Goel for computing the flow in a circular pipe 
having a sudden expansion. Computed values of the Nusselt numbers were 
generally close to the experimentally determined data of Amano et al. 1 

A better modeling of the turbulent diffusive transport processes of 
k and e could be attained by not resorting to the Boussinesq's eddy 
viscosity assumption but instead, by evaluating the diffusive transport 
in terms of the Reynolds stresses in their original form. Thus, the 
modified versions of k- and e -equations can be expressed as 


Dk 
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where C. and C are the diffusion coefficients for k and c, which have 
k e 

been originally recommended to be equal to 0.30 and 0.15, respectively. 
However, upon further investigation, these coefficients have been 
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recommended to be 0.10 and 0.30, respectively in this study. The 
generation rate, G, is accordingly evaluated by using the expression 
given by equation (2.1.12). 


2 . 3 Transport Equations for the Third Moments (Triple Velocity Products) 
of Turbulence 

As discussed in Section 2.1, the turbulent diffusion of the Reynolds 
stresses is governed by the gradients of the third moments (equation 
(2.1.6)). This is expressed, employing Daly and Harlow's 17 model, in 
terms of the product of the Reynolds stresses and their gradients 
(equation (2.2.13)) . 

In order to obtain a transport equation for u^u^u^, t * le transport 
equation for a single fluctuating component of velocity u^ is multiplied 
by u^u^ anc * added to the two other forms of this equation obtained by 
interchanging the indices i, j, k, and finally time averaging as follows 


Du.u. u ^ 
Dt 
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As shown by Chou, 15 this leads to the following transport equation for 
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Equation (2.3.2) represents the transport of the third moments, u.u.u, , 

1 j K 

and shows how their convective transport is balanced by other terms 
which represent the combined transport effects due to generation, 
diffusion, pressure stress and dissipation. 

Terms I and II represent the production rates of u^u^u^ due to the 
action of mean strains and the shear stresses. They need no further 
modification being rather explicit in their form. 

Term III may be rearranged as 


Term III - - (u .u.u, u.) 

9x„ 1 J k 


(2.3.3) 


indicating its diffusive nature. 
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Further simplification of equation (2.3.3) is adapted from Hanjalic 
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and Launder in which the quadruple velocity correlations are expressed 
in terms of Reynolds stresses as 


YjVj ‘ “iWi + u iV u j u i + 'YYVj <2 - 3 - 4 > 

Differentiating equation (2.3.4) with respect to x^ and adding it to 
Term II yields 


3u.u, 3u,u. 3u.u. 

Term II + Term III - -[u u — ^ + u.u — — + ilu. — L ] (2.3.5) 

Equation (2.3.5), hitherto, represents the production rate of u^u^u^ due 
to the interaction of Reynolds stresses with their gradients. Term IV 
indicates the pressure-stress effects and contains the combined effects 
of generation and redistribution. The modeling of the pressure -stress 
term is analogous to the modeling of the pressure -strain redistribution 
term of the Reynolds stresses in equation (2.1.8). 

Term IV is evaluated by first forming the Poisson equation for fluc- 
tuating pressure p by differentiating a transport equation for u^ with 
respect to x^. Upon changing the subscript j to i and the subscript i 
to k yields 
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Integrating this equation with respect to x^ once and then multiplying 
it by u ^ u j a n d finally time averaging the resulting equation yields 
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This integral is approximated as 


d au 

■ • [u i u j 7" <Vj • Vj> + 2 u i u j u i — 1 <2 - 3 ' 8) 

£ dx £ 

Further simplification of equation (2.3.8) is done by assuming it to be 
equal to 

e au 

- C A U i"j\ 7 + C B U i U j"t 7" (2 ' 3 ' 9) 

K 

where and are arbitrary constants . Applying the approximation of 
equation (2.3.9) to the rest of the components in the term IV of 
equation (2.3.2), the final value of Term IV may be expressed as 


Term IV - - u^u^u^ 


(2.3.10) 


where the second term in equation (2.3.9) along with its coefficient is 
merged with Term I in equation (2.3.2) and the value of the coefficient 
C modified accordingly to take up the deficiencies if any in making 
this approximation. 

Term V containing the viscous effects can be rearranged and written 
as 


3 Su.u.u, 3u. 3 cl 

■ Term V - u ( — — ^ — ) - 2i/[u^( — ^ ) 

dX £ 8x £ 1 dX £ dX £ 


3u. 3u, 3u. 3u. 
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J 3x i 3x^ dXj dx £ 


The first part in equation (2.3.11) represents the viscous diffusive 
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dissipation of u u il , 

1 J K 

The modeling of this viscous dissipation term has been adopted from 
the modeling of the dissipation term c ^ in equation (2.1.7). Owing to 
the similarity of the terms in equation (2.3.11) and the dissipation term 
(Term II in equation (2.1.4)) the dissipative effects influencing the 


transport of u^u^u^ can be expressed as 


Dissipation - - (6^. + 8^ + 5^) ck 


1/2 


(2.3.12) 


where the single fluctuating velocities are set proportional to k 


1/2 


with the constant C thrown in to account for the discrepancies, if 
any, in making this correlation. 

Based upon parametric tests, the value of is recommended to be 

0.10. Equation (2.3.2) may then finally be cast in accordance with the 
Reynolds-stress equation (equation (2.1.5)) as 
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(2.3.14) 


where C in P,., , is used to add the second term in equation (2.3.9) 

p 1 j K 

Solution of the transport equations of the third moments (equation 


(2.3.13)) would yield the values of u.u.u, which when plugged into 

l j Jc 

equation (2.1.5) would give the diffusion rate of the Reynolds stresses 
enabling a better approximation and modeling of the Reynolds-stress 
closure. 

2 . 4 Low-Revnolds Number Model for the Third Moments 

In order to incorporate the viscous effects that are predominant 
near the wall, a low-Reynolds number modification has been developed 
here to improve the agreement with experimental data. 


The transport equation for u^u^u^. incorporating this low-Reynolds 
number model is 


D VVk P,„. . F„ 


Dt 


ljk l + ^ k 2 + *ijk + *ijkw ‘ € ijk + D ijk 


(2.4.1) 

This equation is identical to equation (2.3.13) except for the term 
^ijkw t ^ ie additional pressure-stress contribution needed to 


distribute and control the levels of u.u.u, in the viscosity 

i j k 


dominated 


near-wall region. 



The theory of this low-Reynolds number model is based upon the 
existence of three regions: i) viscous sublayer, ii) buffer layer, and 


iii) the fully turbulent core. The viscous sublayer exists upto a 
distance y + < 10.0 where 


+ f> <V p) y 

y 
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In this region the dissipation rate is governed by the equation 
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The dissipation rate in the buffer layer is formulated as 


dissipation - C 
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Thus, the terms + ^ijkw are to S et ^ er given as 
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where C and C are modified coefficients for the pressure -stress and 
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its corresponding near-wall approximation. The symbol | A. , B | 
represents the largest of the arguments in this bracket. 

The development and incorporation of equations (2.4.2) and (2.4.3) 
is based upon the near-wall model of Amano et al ., 1 where boundary 
conditions of the wall have been derived using these profiles for the 
dissipation rate. 

In equation (2.4.3), is another constant defined as 

C. - *C _3/4 (2.4.5) 

** £ \i 

where k is the von Karman constant and is equal to 0.42. 



Figure 2.1 shows the variation of the dissipation rates obtained by 
equations (2.4.2) and (2.4.3), with the non-dimensional distance y + from 
the wall. The dissipation rate, e, obtained by the solution of the e- 
equation (equation(2 . 2 . 6) ) is also shown for comparison. 

As shown in the figure, the dissipation rate is considerably 
augmented in the near-wall region and therefore justifies the inclusion 
of equations (2.4.2) and (2.4.3) for correcting the near-wall values of 
U i U j U k t * 16 low -Reynolds number model. 

The dissipation rate given by equation (2.4.2) is higher in the 
near-wall region but decays very rapidly with increasing distance from 
the wall, as seen in the figure. Thus, the incorporation of equation 
(2.4.4) in the transport equations of the third moments enables the 
viscosity dominated near-wall region to exercise its influence on the 
behavior of the third moments appropriately. 

2.5 Algebraic Equations for the Third Moments of Turbulence 

Algebraic correlations have been developed by different researchers 
to evaluate the third moments, without solving the complete transport 
equation for u^u^u^. 

In order to arrive at these algebraic equations considerable 
simplifying assumptions have been made which tend to introduce imperfec- 
tions in the modeling. 

Based upon their experimental measurements of asymmetric flow in a 

25 

plane channel, Hanjalic and Launder found that the term i n 

equation (2.3.13) was negligible compared to other terms in the 
equation, especially for wall boundary layers. Furthermore, upon 
neglecting the convective transport as well as the viscous dissipation 
and molecular diffusion, the equation (2.3.13) reduces to 
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(2.5.1) 


Upon substituting the values of ^ and P._. v from equation (2.3.14) 


' ijk. 


into equation (2.5.1), Hanjalic and Launder obtained the following 
equation 
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Daly and Harlow 17 simplified equation (2.5.2) by suitable contraction 
and obtained the following correlation 


(2.5.2) 
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After further simplification, Shir 46 proposed 


u.u.u. --0.04 =— 1- 

i J Tc 


k du^u. 
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(2.5.4) 
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Cormack et al. obtained these algebraic expressions by approximating 


the experimentally determined profiles for u.u., e and u u u. with 

1 J i j K 

polynomials with coefficients, chosen to give a least squares fit to the 

data and discretized the cross -stream coordinate in each flow into a 

certain number of grid intervals. They used the most general model for 

the triple velocity correlation tensor as generated using the asymptotic 
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approach of Lumley and Khajeh Nouri. 

They obtained a set of algebraic equations and used the least 
squares approach to solve these set of equations and optimized the 
parameters. The model they obtained is as follows 
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Out of the 20 parameters (a^) that they had started out with, they 
were able to determine the most significant ones as a^, a^, a^, and a^. 
The values of these parameters have been recommended for various kinds 
of flows along with the universal values applicable to all flows wherein 

e x - *8.14 x 10' 3 
o 5 - -1.72 x 10* 2 
a 7 - -4.80 x 10‘ 2 

a n - -1.02 x 10' 1 (2.5.7) 


In arriving at the above mentioned algebraic correlations (equations 
(2.5.2) - (2.5.5)) various simplifying assumptions had been made which 
necessitated the deletion of different terms from equation (2.3.13) 



thus increasing the deviation from the actual physical phenomenon. A 
more realistic modeling would be obtained if the complete transport 
equation of the third moments could be solved iteratively. 

Table 2.1 lists the recommended values for constants used in the 
present study. 


Table 2 . 1 Recommended Values for the Constants Used In Turbulence 
Modeling . 
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B 1 

- (C 2 + 8)/ll 

B 2 

- (30C 2 - 2)/5 

B 3 

- (8C 2 - 2)/ll 

C 1 

- 1.5 
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CHAPTER 3 


NUMERICAL PROCEDURE 


3 . 1 Formulation and Discretization of the Transport Equations 

The steady state 2-dimensional equations for all the parameters 
solved in this study can be, in general, cast and expressed as 


8pV<f> 8pV<f> 8 d<f> d 8<f> 

+ [— (r. — ) + — (r — ) ] + s - o 

ax ay ax x ax ay * ay 9 (3.1.1) 


where «j> stands for different dependent variables U, V, k, e, u ^ u j and 
u^u^u^. r x and ^ are the corresponding diffusion coefficients in the 
streamwise and the normal directions to the flow respectively. The re- 
maining terms are lumped into the source (or sink) term S,. 

9 

Equation (3.1.1) shows how the convective transport (net rate of 
inflow of <f> to the cell by convective fluxes) is balanced by the 
diffusive transport (net rate of inflow of <f> to cell by diffusive 
fluxes) and the generation or dissipation rate of <f> within the cell. 
Discretization of equation (3.1.1) is carried out by dividing the 
computational domain into elemental control volumes and then integrating 
the equation over each of these control volumes. Figure 3.1 shows a 
typical node P inside its control volume bordered by neighboring nodes 
E, W, N, S (representing the east, west, north and south locations). 

AO 

Following Patankar, equation (3.1.1) is integrated over this 
control volume and thus can be written as 


Vp " Ve + Vw + Vn + Vs + s u 


(3.1.2) 


where 


a E * V ( I P .I ) + M- F ..°M 

<Hj - V ( l P wl> + < 3 ' X ' 3 > 

*S- V ( l P nl> + ll-F n .O|J 

a s ‘ V ( l p « l> + tl F ,.°l) 

“P ‘ a E + “B + a K + a S • S P <3 ' 1 ' 4> 

and 

JS^Vol - + S p * p 

In equation (3.1.3) A(|P^|) is a function of the cell Peclet number P^ 
whose formulation depends upon the type of differencing scheme being 
employed. The aspect of differencing scheme is discussed in the 
following section. The symbol [|a,B|] is the same as the one used in 
(2.4.4). 


Equation (3.1.2) is the discretized algebraic form of the original 
differential equation (3.1.1). 

Figure 3.2 shows the staggered grid system employed in this study. 

The scalar cell denoted by the node P is used to discretize and compute 

2 2 " 1 

the scalar parameters like P, k, e, u , and v and u^u.u^. 

The U-momentum cell is staggered with respect to the scalar cell by 
shifting it by half a cell to the left while the V-momentum cell is 
staggered by shifting it half a cell below the scalar cell. The cell 
for uv, however, is staggered by shifting it half a cell to the left and 
half a cell below the scalar cell as shown in the figure. 

The reason for doing this staggering of the cells is the advantage 
gained by making the velocities available at locations where they would 
facilitate gradient calculation thus making it easier to satisfy the 
continuity equation. Likewise, the stresses are evaluated at locations 
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where they would enable their gradients and production rates to be 
evaluated easily. 

3 .2 Differencing Scheme for the Discretized Equations 

After having transformed the original differential transport 
equation (equation (3.1.1)) into an algebraic equation (equation 
(3,1.2)), it is required to prescribe a formulation for the combination 
of the convective and diffusive fluxes such that the coefficients in 
equation (3.1.2) may be evaluated. Formulation of a scheme to combine 
the convective and diffusive fluxes is referred to as differencing 
scheme. 

The differencing scheme employed in this study has been adopted from 
Amano et al. 1 in which the exponential (exact) scheme in Table 3.1 is 
expanded upto the fourth order term as 

1 1 0 1 * 

a(|p i) - [ i o , (i - - |p i + — |p r - — |p .ni] (3.2.D 

1 2 1 12 720 

As shown by Amano el al . , 1 the deviation of the hybrid scheme from 
the exact exponential scheme is quite large at P^ - ± 2 and so it is 


TABLE 3.1. Function A( PjJ) for Different Differencing Schemes 


Scheme 


Function A( P^) 


Central differencing 

Upwind 

Hybrid 

Exponential (Exact) 


1.0 - 0.5 |P | 

1.0 

[ 1 0 , (1 - 0.5 |P.|)|] 

|p|/(e l P iI - 1) 

1 1 

1 1 o , (i - -IpJ + — |PJ 

2 1 12 1 


720 


pj 4 ) 


Fourth Order 
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rather premature to neglect the diffusion effects as soon as |P.J 
exceeds 2 . 

The fourth order scheme of equation (3.2.1) reduces to the upwind 
scheme and neglects the diffusion effects after |P.J becomes greater 
than 4 which makes it approach closer to the exponential (exact) 
solution. Moreover, the fourth order scheme is less expensive to 
compute compared to the exponential scheme, although it may involve 
extra expenses compared to computation with the hybrid scheme. An 
additional advantage of the fourth order scheme lies in the fact that 
unlike the hybrid scheme, it introduces very little or no false 
diffusion effects into the solution. False diffusion is an inherent 
aspect of the hybrid differencing scheme making it less accurate at high 
Reynolds number flows, convection dominated flows and more so, for 
recirculating flows. 

3.3 Solution of the Equations 

After having formulated and descretized the equation as shown in the 
previous sections, they are solved iteratively using line-by-line 
solution method. The initially guessed values of the flow field for 
various parameters is prescribed and these are improved upon from one 
line to the other in successive iterations. 

M2 

As shown in Roache, for the solution of the equations along the 
nodes on each line (e.g. North- South line), values of the parameters on 
the neighboring lines are assumed to be temporarily known (previously 
stored values). The equation for each node on this North-South line 
then reduces to one where only three values (^ p , ^ and ^ g ) are unknown. 
The set of equations for all nodes, on the North-South line then takes a 
simple form in which the non-zero coefficient matrix is tri- diagonal and 
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which is solved by the line-by-line iterative technique called TDMA 
(Tri-Diagonal Matrix Algorithm). 


3 . 4 Boundary Conditions 

The numerical solution of the transport equations requires the 
provision of correct boundary values in the flow domain. 

For backward- facing step flows, boundary conditions need to be 
prescribed at the walls and along the outflow section at the exit. 

At the wall boundaries, the "wall law" is used to specify mean 
velocities, turbulence kinetic energy, the dissipation rate, the 
Reynolds stresses and the third moments. 

For the mean velocities, the wall shear stress is determined using 
the velocity gradient at the wall, which is then used to prescribe the 
boundary values . 

For the turbulent kinetic energy, the generation rate at the wall is 
modified by incorporating the wall shear stress in the generation rate, 
which is then substituted in the transport equation of k thereby 
introducing the effect of the wall. 

The dissipation rate, however, is expressed in terms of the kinetic 
energy as 


3/2 




(3.4.1) 


which is similar to the equation (2.4.3). This value constitutes the 
boundary value for the dissipation rate. 

For prescribing the wall values for the Reynolds stresses, Launder 


et al. 3A had obtained proportionality correlations between u^Uj and u r f° r 


channel flows as 
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2 y l dP 

u.u. - C..U - (1 - 6..) — 

i j li r iv , 

J J p dy 


(3.4.2) 


where U is the friction velocity and where y. and y_ represent the 

T 1 Z 

coordinates normal and parallel, to the wall in consideration, 

respectively. The coefficients C^. are given as 


c n - 5 - x 
C 12 “ - 1 ' 0 
C 22 " 1,0 


(3.4.3) 


Using the correlation between k and in the wall proximity region, 
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derived by Hanjalic and Launder, as 


k « 3.5U 


(3.4.4) 


the boundary values of equation (3.4.2) may be expressed as 


y-. dP 

u.u. - C..k - (1 - 8 ..) 

l j ll lj , 

J J J P dy 2 

where the coefficients are accordingly changed as 

c u “ 1 - 214 

C 12 - -0.238 
C 22 - 0.238 


(3.4.5) 


(3.4.6) 


For prescribing the wall values for the third moments, the algebraic 
correlation of Shir/ 8 equation (2.5.4), combined with the wall values of 


u^Uj , equation (3.4.5), is used as 


k 2 d 




y x dp 


- 0.04 {C..k - (1 - 6..) — ) (3.4.7) 

« ^ ^ P dy 2 


It was seen that the solution of the equation for the third moments did 


not seem to be affected by the choice of boundary conditions; for ins- 
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tance, even a value of zero at the wall did not change the solution 
appreciably. 

For prescribing the boundary values along the outflow section, at 
the exit, a condition of continuous flow with zero streamwise gradients 
was imposed for all the parameters. 

3 . 5 Computing Details 
a) Grid Selection 

In order to have a fine mesh in the recirculation and near -wall 
regions of the computational domain, a variable grid system was employed 
wherein the grid spacing increased monotonically in the downstream flow 
direction thus generating a fine spacing near the step and a relatively 
coarser mesh towards the far downstream end. In the transverse 
direction, however, the grid spacing was made finer near the bottom wall 
and it became coarser with increasing distance from this wall. 

Exploratory tests were made using mesh sizes of 32 x 32, 42 x 42 and 
52 x 52 with different grid expansion factors in the streamwise and 
transverse direction and with different axial lengths of the solution 
domain. Amano and Goel 5 and Sindir 50 employed the change in 
reattachment length as a basis for determining the grid independency 
stage and thus came up with a 42 x 42 mesh size to be grid independent. 

The study carried out here revealed that the shift in the reattachment 
length with changing the grid size and spacings is small compared to the 
change in other parameters like the mean velocities and the Reynolds 
stresses at different regions of the domain, making it difficult to 
ascertain a criterion for grid independency tests. Tahry 53 used the 
Reynolds -Stress Model with a 32 x 30 grid to predict flow in an engine 
cylinder and observed a shift in the mean velocity profiles at 144°ATDC 
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when a 45 x 45 grid was employed. According to Armaly et al. 6 the 
required grid number to ensure grid independence and the number of 
iterations (for the same grid number) increases with the Reynolds 
number. They could not attain complete grid independence with a grid 
number of 2600 (51 x 51 grid). Likewise Chieng observed a marked 
improvement in the Nusselt number predictions for flow into a pipe with 
a sudden expansion when the grid was changed from 22 x 22 to 50 x 50. 

In order to determine a reliable grid system, the grid size was kept 
fixed at 52 x 52 and exploratory tests were carried out with different 
grid spacings and lengths of the solution domain. Streamwise mean 
velocity and Reynolds stress profiles were compared with the 
experimental data for that particular set of experimental flow 
conditions. The combination of the grid parameters which had the best 
agreement was selected to be the final grid system for that particular 
set of experimental flow conditions . Specific grid details will be 
mentioned in Chapter 4 while comparing the results of the computation 
with experimental data. 

b) Convergence Criterion . 

After having specified the grid system, the cell dimensions and 
distances are evaluated in the computer program. Initial values of the 
parameters are then prescribed prior to starting the iteration process. 

The convergence criterion of the solution is based on the conserva- 
tion of mass and momentum in the flow domain. Formulation of equation 
(3.1.2) for the mean streamwise and transverse velocity components 
expresses the momentum balance in the flow domain. Satisfaction of 
these equations upto 1% of the total momentum inflow at the inlet of the 
channel, constitutes the criterion for convergence. Likewise, 
formulation of equation (3.1.2) for pressure, expresses the mass balance 
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in the flow domain. Satisfaction of this equation upto 1% of the total 

mass inflow at the inlet of the channel, also constitutes the criterion 

for convergence. The total absolute residual source values for each of 

these equations is monitored for every iteration and the satisfaction of 

the above three criterion is taken to be the completion of convergence. 

In the solution of the transport equations of the third moments, the 

total residual source terms for each equation are normalized by 

across the inlet, and are monitored for each iteration. Convergence of 

the solution for u^u^u^ is assumed to be completed if either the 

-12 

normalized maximum total residual source is less than 5.0 x 10 or if 

the maximum difference between the residual source terms of two 

-12 

consecutive iterations is less than 3.0 x 10 

In the solution of equations of the RSM incorporating the third 
moments, the above mentioned convergence criterion needs to be satisfied 
to qualify for the completion of convergence. However, the 
consolidation of the RSM by incorporating the transport equations of the 
third moments requires a large amount of computations . In order to 
conserve computational time , the convergence criteria for the third 
moments is made less stringent by relaxing the above mentioned limits to 
1.0 x 10~’ L ° and 1.0 x 10 respectively, 
c) Model Sequencing 

The strong nonlinearity of the system of the equations makes 
this solution sometimes unstable and hence, susceptible to divergence. 

In order to avoid divergence, the new models developed in the study are 
introduced into the solution in sequence. 

The solution is begun with the standard k-e model. When the absolute 
residual source values drop down to about 3%, the Reynolds -Stress Model 
is invoked which entails the solution of the transport equations of u ^ u j 
along with the momentum and k-e equations iteratively. 
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In the next sequence, the modified k and e equations are exchanged with 
the original version when the source term is down to about 2%. The 
third order closure is introduced in the final sequence when the source 
term is 1.5%. Convergence is assumed to be complete when the absolute 
residual source term drops down to 1% and when the source terms for 
u.u.u. have met their convergence criterion. By sequencing the solution 

1 j K 

of the equations (models) the solution process became less prone to 
divergence . 

Specific grid and model sequencing details are mentioned in Chapter 
4 while comparing the computational results with experimental data. The 
computations were carried out using a UNIVAC 1100 computer. The CPU 
time and iteration number/history details vary with the turbulence 
models and with the input data, and are hence mentioned in Chapter 4 
along with the specific grid and model sequencing details. 
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CHAPTER 4 

RESULTS AND DISCUSSION OF COMPUTATIONS IN BACKWARD -FACING STEP FLOWS 
4 . 1 Introduction 

This chapter presents the results of the study performed in 
mathematically modeling the turbulence phenomena of separating and 
reattaching flows. The models developed and presented in Chapter 2 are 
applied to plane backward- facing step flows and are compared with 
existing experimental data. 

As mentioned earlier, the reason for selecting backward- facing step 
flows to test the models is that this class of flows possesses most of 
the important features of flow separation and reattachment. They are 
greatly susceptible to modeling modifications but, on account of their 
simple geometry, they are easier to program and compute. Moreover, 
there exist considerably large sets of reliable and recent data on this 
class of flows which facilitates better comparison with the modifica- 
tions in modeling. 

In the previous paper 5 , the results obtained with the RSM by solving 



the results obtained by the k-c model, in which the Boussinesq 
correlation, equation (2.1.1), is used to evaluate the Reynolds 
stresses. Comparison was also made with the computed results of 

50 § 

Sindir using the Algebraic Stress Model (ASM) and the experimental 




The ASM was developed by Rodi by relating the convective and 
diffusive transport of the Reynolds stresses to the turbulent kinetic 
energy k, which after simplification and rearrange ment reduce to a set 
of algebraic equations for evaluating the stresses u^u. . 
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data of Seegmiller and Driver and Kim et al . who had employed a 
laser doppler velocimeter (LDV) and hot wire anemometer, respectively, 
to acquire their data. 

Regarding agreement with the experimental data, the results obtained 
by the RSM showed the best agreement in the recirculating region for 
predicting normal stresses and showed the best agreement in the 
reattachment region for predicting shear stresses. 

In general, it was found that the k-e model underpredicts whereas 
the ASM tends to overpredict the levels of the shear and normal 
stresses. The results of RSM, however, tend to remain in between the 
results obtained by the k-e model and the ASM. 


4.2 Transport Equations for Third Moments of Turbulence with Low- 
Reynolds Number Model 


Third moments evaluated by solving the transport equations of u.u.u^ 

are compared with the experimental data of Driver and Seegmiller 1 

12 


19,20 


and Chandrsuda and Bradshaw. 

After obtaining a good agreement with the experimental data in 
consideration, the final converged values of the mean velocities, 
Reynolds stresses, k and e are stored and are used later to solve the 


transport equation for u^u^u^ iteratively. 

In the high-Reynolds number model of the transport equation for 

u^UjU^, the value of the empirical constant for the pressure -stress 

effects e Q uat ^ on (2.3.10) is recommended to be 5.8, after 

matching the peak values of the third moments with experimental data. 

In the low-Reynolds number model, the values of the empirical 

constants C and C for the combined effects of the pressure -stress 

7 yw 
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effect and the near-wall effect, 0.^ + ^ijkw’ were su i ta kly modified 
and were set equal to 2.0 and 8.0, respectively. 

However, in solving the transport equations by both the high-as well 
as the low-Reynolds number models the effects due to the dissipation 


rate, e. , defined by equation (2.3.12) have not been incorporated for 
ij Jc 

the results presented in this section. This is because the results 


presented in this section are for reference only. The final results 


presented in the subsequent section have been obtained by the 


consolidation of the RSM and contain all the effects discussed earlier. 


Figures 4.1 and 4.2 show the computed results of the third moments 
obtained by both high- and low-Reynolds number models of the transport 
equation and compare them with the two sets of experimental data 
mentioned earlier. It is seen that both the high and low-Reynolds 
number models give similar levels except near the wall, wherein the 
results obtained by the low-Reynolds number model are improved and are 
considerably superior. As discussed earlier, the reason for this 
improvement is attributed to the enhancement of the diffusive transport 
in the near-wall viscous region due to the additional dissipative 
effects . 


Figures 4.3 and 4.4 compare the profiles of the third moments 
obtained by the four algebraic models (equations (2.5.2), (2.5.3), 
(2.5.4), and (2.5.5)) and by those obtained through the solution of the 
transport equations, incorporating the low-Reynolds number model, with 
the experimental data. As shown in these figures, the results by the 
transport equations for the triple velocity correlation give a much 
better agreement with the experimental data than the results obtained by 
the corresponding algebraic models. The transport model also improves 
the prediction near the wall. 


43 


ORIGINAL PAGE 9 

The transport model is very versatile in its 

application; its predictions are reliable in both homogeneous and 
inhomogeneous flows on account of its symmetry property in all the three 
directions . 

The values of the emperical constants , C and C were modified 

7 7W 

slightly after the incorporation of the dissipation term e.., and after 

ijk; 

merging the transport equation model for u.u.u, with the RSM. Details 

l j R 

of this follow in the next section. 


4.3 Consolidation of RSM 


Results of the third moments presented so far had been obtained 
implicitly by retrieving and using the converged and stored values of 



without resolving the equations of those other parameters. 


In the consolidation of the Reynolds -Stress Model, however, the 
transport equation for all the parameters like U and V velocities, k, e, 
u^u. and the third moments, u^u^u^, are solved together iteratively. 

This system of equations thus involves the iterative solution of a total 
of 11 equations all of which are nonlinear and intercoupled. 

As mentioned in the preceding section (4.2), the values of the 
empirical constants and had to be suitajj^r modified during the 
consolidation of RSM in order to comply with the new system of 
equations. The final recommended values of these constants are 3.0 and 
8.0, respectively . 

The solution is initially begun by solving the momentum equations 
with the original k and € equations (equation (2.2.5) and (2.2.6)). 

Model sequencing is employed to introduce the new equations into the 


solution stream. 
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When the absolute residual source term reaches 3%, the transport 
equations for the Reynolds stresses are introduced into the computation 
stream and the iterative process is continued. At the source term value 
of 2%, the new equations for k and e with the modified diffusive 
coefficients (equations (2.2.7) and (2.2.8)) are introduced and 
iterative process is carried on. At the residual source term value of 


1.5%, the transport equations for u.u.u, , (equation (2.4.1)) are finally 

1 J R 

merged into the computation stream. Complete convergence is assumed 
when the absolute residual source term drops to less than 1%. 

19 20 

Experimental data of Driver and Seegmiller ' and Chandrsuda and 
12 

Bradshaw are compared with computational results obtained by three 
turbulence models explained hereunder. 

Model I is the standard k-6 model where the Reynolds stresses are 
evaluated by using the Boussinesq's approximation. A total of 4 
transport equations get solved in this case. 

Model II refers to modified k-c model that the Reynolds stresses are 


obtained by solving the transport equations for u ^ u j which contain the 


17 


algebraic correlation of Daly-Harlow for u^u^u^. A total of 7 
transport equations get solved here. 

Model III, however, uses the modified equations for k and e, and 


solves the transport equations for u ^ u j as well as u^u^u^. A total of 
11 transport equations are solved iteratively in this case. 

In the following figures, unless specified otherwise, the following 
legend is to be used to identify the models. 

Model Line Type 

Model I 

Model II 

Model III 


O 


Experimental Data 
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Figures 4.5 and 4.6 compare the U velocity profiles obtained by the 
different models with the two sets of experimental data. 

As shown in figure 4.5, the results obtained by Model III are 
generally closer to the experimental data. The agreement improves in 
the shear layer and is more encouraging in the redeveloping region. The 
three models generally give similar results though the agreement is not 
quite as good with the data in figure 4.6. 

The similarity of the results obtained by the three models is 
attributed to the fact, that in the solution of the momentum equations, 
the Reynolds-stress terms are being handled through the turbulent 
viscosity approximation, and so do not appear in those equations 
explicitly. This is because, due to the greater significance of the 
momentum equations (mass balance for proper convergence) , they tend to 
be very sensitive to modeling modifications thereby making the solution 
unstable. In order to overcome this instability and to avoid 
divergence, the turbulent viscosity (Boussinesq' s approximation) is 
retained in the momentum equations . 

Figures 4.7 and 4.8 compare the shear stress, uv, profiles. It is 
seen that for any given axial location, all the three models predict a 
more or less similar transverse location of the peak values. The levels 
predicted by Model I are the lowest of all the three models . Model II 
tends to overpredict in figure 4.7 whereas Model III overpredicts in 
figure 4.8. It is seen that the predictions by Model III have generally 
a better agreement with the experimental data. 

2 

Figures 4.9 and 4.10 compare the normal stress u profiles at 
different axial distances from the step. The predictions by Model I are 
always lowest whereas Model III predicts the highest values. 
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Predictions by Model II seem to give a better agreement with the 
experimental data. 

Figures 4.11 and 4.12 compare the normal stress v profiles at 

different axial distances from the step. It is seen that Model I 

exhibits an opposite trend when predicting the normal stresses u and 
2 

v . Model I predicts the largest values whereas Model II predicts the 
lowest values. The agreement with the experimental data in and around 
the recirculation region is better with Model III. 

Figure 4.13 compares the profiles of the kinetic energy k obtained 
by different models. Experimental values of the kinetic energy have 
been obtained by multiplying the sum of the normal stresses (u and v ) 
by 0.75. It is seen that Model III gives the best predictions with 
the experimental data whereas Model I is quite deficient in its 
predictions . The improvement made in the RSM due to the incorporation 
of the transport equations of the third moments and the modified 
equations of k and e seems to be encouraging from this figure. 

Figure 4.14 compares the energy dissipation rate profiles. Model I 
seems to give a better agreement whereas Model II is not quite as good. 
Predictions by Model III are intermediate. 

Success in predicting the trends of the various parameters discussed 
above not only depends upon the type of model being employed but also 
upon the particular parameter being considered and the location inside 
the flow domain. 

The choice of models had very little influence on predicting the 

mean velocity profiles. In predicting the Reynolds stresses, however, 

Model III gave superior results in and around the reattachment region. 

2 

Model I seemed to give better predictions for v and energy dissipation 
rate profiles. Results obtained by Model II seemed to be closer to the 
results by Model III and generally were better than those by Model I. 



47 


Figures 4.15 through 4.21 show the profiles of the various 
components of the third moments at different axial locations from the 
step obtained by Models II and III. 

Profiles obtained by Model III are consistently superior to those 
obtained by Model II. The discrepancy with the experimental data is 
rather large for the profiles of vw; both the models underpredict the 
levels considerably although Model III is still better than Model II. 

As discussed earlier, the solution of the complete transport 
equations of the third moments entails taking into consideration, the 
phenomena of convection, diffusion, generation, redistribution and 
dissipation which consequently yields greater accuracy in the turbulence 


model . 
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Table 4.1: Computational Details 


Computational Details 


MODEL EQUATIONS Driver and Seegmiller 19 Chandrsuda and Bradshaw 12 



BEING 

SOLVED 

(Yo/H - 8 

.0) 

(Yo/H - 2. 

• 5) 


Iteration 

Number 

CPU 

Time 

(mins) 

Iteration 

Number 

CPU 

Time 

(mins) 

I 

2.1.2 

2.2.5 

2.2.6 

265 

25 

500 

36 


2.1.2 





II 

2.2.5 

2.2.6 
2.1.15 

222 

28 

443 

42 


2 . 1.2 

2.2.7 

III 2.2.8 240 46 

2.1.15 
2.4.1 


603 


72 


CHAPTER 5 


CONCLUDING REMARKS 


The salient feature of this study lies in the fact that it has taken 
the Reynolds Stress Model of turbulence to a higher degree of 
sophistication and accuracy through the development and subsequent 
introduction of the transport equations of the third moments into it. 
Being primary diffusive transporters, the correct simulation of the 
third moments, thus enables genuine government of the Reynolds stresses 
in the flow fields with separation, reattachment and recirculation. 
Inadequacy of the existing algebraic equations for evaluating the third 
moments was demonstrated, and approximate correction factors were 
recommended for these equations by comparing them with the experimental 
data. 

Development of the transport equations of third moments required the 
formulation of their equation and then the modeling of the various terms 
in the formulated equation which was sometimes attained through the use 
of empirical constants. The relative effect of the constants on the 
behavior of the third moments was investigated and suitable values of 
these constants were recommended. 

A low-Reynolds number model was developed and incorporated into the 
transport equations of the third moments which greatly improved their 
predictions near the wall. 

With regard to the Reynolds stresses, however, the RSM predicted 
better normal stresses in the recirculating regions and better shear 
stresses in the reattachment region. 



The highlight of this theoretical study is the consolidation of the 
RSM by merging it with equations of k and e containing the new 
recommended diffusion coefficients and the transport equations of the 
third moments containing the low-Reynolds number effect. This system of 
equations required the iterative solution of 11 non-linear and 
intercoupled equations. The mean velocity, Reynolds stresses, k, e and 
third moments profiles thus obtained were compared with two sets of 
experimental data and with similar computational results obtained by the 
k-e model and the original RSM. 

In the modeling process and the determination of empirical constants 
the experimental data have been relied on heavily. Since the experi- 
mental data are . associated with various levels of uncertainties, the 
corresponding theoretical models are also associated themselves with 
some of these uncertainties which may be reduced in future studies. 
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